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SUMMARY 


The effecte of power on the stability and the control 
eharacteri sties of an airplane are discussed and methods 
of analysis are given for evaluating certain dynamic char- 
acteristics of the airplane that are not directly discern- 
ible from wind-tunnel tests alone. Data are presented to 
show how the characteristics of a model tested in a wind 
tunnel are affected hy power. 

The response of an airplane to a rolling and a yawing 
disturbance is discussed, particularly in regard to changes 
in wing dihedral and fin area. Solutions of the lateral 
equations of motion are given in a form suitable for direct 
computations. An approximate formula is developed that 
permits the rapid estimation of the accelerations produced 
during pull-up maneuvers involving abrupt elevator deflec- 
tions. 


INTRODUCTION 


Some time ago, the NACA undertook an investigation to 
determine the flying qualities of a low-wing, pursuit mon- 
oplane. This airplane (fig. 1), like many of its type, 
was found to possess several undesirable flying character- 
istics. Accordingly, a number of modifications to the 
airplane were recommended. In order to study the effects 
of the proposed modifications and in order to investigate 
the reliability of wind-tunnel data for estimating the 
behavior of the full-scale airplane, a scale model of the 
airplane, with and without the proposed modifications, 
was tested in a wind tunnel. Tests were made with and 
without the propeller operating in order to study the ef- 
fect of power on the characteristics of the model. 
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During the course of these wind-tunnel tests, sev^rM 
Interesting results were obtained concerning the effects <> 
power on the characteristics exhibited by the model. In 
the evaluation of the data, certain methods of analysis 
were employed that were found useful In estimating the 
dynamic behavior of the airplane, particularly as regards 
the effects of the proposed modifications. These subject* 
are treated in the present paper. Zt is felt that this 
paper will serve as an aid to other wind-tunnel investiga- 
tors seeking to evaluate the flight characteristics of an 
airplane from wind-tunnel results. 

The first portion of the paper is concerned with ths 
effect of power on certain of the wind-tunnel character- 
istics exhibited by a model. Graphs are presented that 
emphasise the importance of obtaining wind-tunnel data 
with powered models. The second portion of the paper d U- 
cusses the lateral motions of the airplane. Concrete a] - 
plications are included showing how the solutions to the 
lateral equations of motion were used to estimate the 
change In dynamic flight characteristics likely to result 
from the use of the recommended modifications to the air- 
plane. finally, the longitudinal motion of the airplane 
is considered. An approximate formula Is developed that 
was found useful in estimating the accelerations and the 
stick forces likely to develop in an airplane during an 
abrupt pull-up maneuver. 


TH1 IMPORTANCE Of POWERED MODELS IN WIND-TUNNEL TESTS 


Although the Influence of the slipstream on the char- 
acteristics of an airplane has been appreciated for many 
years, wind-tunnel tests with propellers operating have 
not been the usual procedure. The frequency with which 
undesirable flying characteristics appear in modem air- 
planes has led to some apprehension concerning the relia- 
bility of wind-tunnel data obtained without powered model* 
The particular teeting technique employed in the present 
investigation will be reserved for detailed discussion in 
a future paper. It will suffice to say here that the 
power-on test procedure should be such that data may be 
obtained for thrust conditions of the model corresponding 
to those of the full-scale airplane at the lift coefficlan 
and the flight attitude under consideration. 
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Figure 2 ie a plot of the elevator angle required for 
trim against the indicated air speed at which trim occurs. 
These curves were evaluated from data secured from (l) 
wind-tunnel tests with no propeller on the model, ( 2 ) ac- 
tual flight tests with the full-scale airplane, and (3) 
wind-tunnel tests with a powered model. The shapes of 
the curves are indices of the static stability (see later 
disousaion of longitudinal motion) displayed by the air- 
plane in flight and by the models in the wind tunnel. The 
agreement between the flight results and the results ob- 
tained with a powered model is good , part icularly for the 
flap-up condition. The larger discrepancy existing for 
the flap-down condition is partly attributable to the 
rather large tab deflection (11° nose up) used on the 
full-scale airplane for this test. The related variation 
of the elevator stick force required for trim with indi- 
cated air speed is not presented here because the model 
data and the flight data were not directly comparable 
owing to the absence of a trim tab on the model. In the 
evaluation of the information obtained from wind-tunnel 
tests, certain methods of presenting the data have been 
found to facilitate the analysis. The resulting pitching 
moment acting on an airplane trimmed for steady equilib- 
rium flight is zero. This condition is true whether or 
not the airplane is flying yawed. As far as the pilot is 
concerned, the dihedral effect of the wing manifests it- 
self as a pure rolling moment about an axis lying in the 
plane of symmetry of the airplane. In wind-tunnel in- 
vestigations, on the other hand, the rolling moment is 
usually measured about an axis coincident with the tunnel 
axis. Consequently, when the model is yawed relative to 
the tunnel axis, the measured rolling moment about this 
axis will contain a component of the pitching moment act- 
ing on the model unless the model is trimmed for zero 
pitching moment. In the comparison of the slopes of 
different rolling-moment curves obtained in yaw tests, it 
is therefore necessary to know whether the measurements 
were made with the model trimmed in pitch. 

If the system of axes discussed in the subsequent 
section on lateral motion is employed for the presentation 
of wind-tunnel data, the necessity for continually check- 
ing the trim of the model during yaw tests is avoided. 

In addition, the data in this form can be used directly 
in the evaluation of the stability derivatives necessary 
for the solutions of the equations of motion. 
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In figure 3 are plotted for various power conditions 
the rates of change with angle of yaw of the yawing-moment 
coefficient dC n /d\|/, the rolling-moment coefficient 
dC^/dty, and the lateral- force coefficient dCy/d ii. The 
decrease in dihedral effect dCj/dty with the application 

of power should be observed. In the condition presented 
(flap up), the decrease in dihedral effect with the applie*- 
tion of power is not great, but, in the flap-down condition 
and at high lift coefficients, the effect of power may he 
critical. 

An inspection of figure 4 will suffice to illustrate 
the effect of the slipstream on the characteristics of a 
control surface placed within its boundaries. The nega- 
tive yawing moment that exists with neutral rudder when 
the propeller is operating should be observed. The moment 
is apparently caused by the rotation of the slipstream. 

It is, however, recognised that the slipstream rotational 
characteristics exhibited by the model in this test are 
probably not identical with the slipstream rotational 
characteristics that would be obtained in flight with the 
full-scale airplane, even under similar thrust condition*. 
The results do indicate, nevertheless, that with the pres- 
ent trend toward greater and greater power, the effect of 
tnis slipstream rotation may be a critical factor in tail 
design. 

Figures 2, 3, and 4 have been included to illustrate 
the effects that power may have on the stability character- 
istics of as airplane and on the effectiveness of the hor- 
isontal and the vertical tail surfaces. These effects are 
difficult te estimate and consequently recourse must be 
had to wind-tunnel tests of a model equipped with running 
propellers. Wind-tunnel tests of a powered model require 
more care and are more expensive to perform than tests of 
the conventional type of model. By a judicious choice of 
tests, however, it is possible to secure the information 
greatest value with a minimum of time and expenditure. 


thi lateral motions of the airplane 


Certain lat eral- stability and control characteristics 
of an airplane can be directly determined from wind-tunnel 
tests of a model. These characteristics usually Include 
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(1) the directional stability (fig. 3(a)), (2) the dihe- 
dral effect (fig. 3(b)), and (3) the static effectiveness 
of the control surfaces (fig. 4). This information is in- 
dispensable to the designer if he is to proportion his 
airplane properly and is also necessary for the evaluation 
of the stability derivatives from which the dynamic be- 
havior of the airplane may be predicted. 

The dynamic flight characteristics that an airplane 
will display are not, however, readily obtained. The be- 
havior of an airplane subjected to a lateral disturbance 
is complicated by the coupling of the ensuing rolling and 
yawing motions. The coupling of these two motions makes 
it impossible to estimate the effectiveness of the ailer- 
ons and the rudder in flight from static tests of the con- 
trols alone, a sound evaluation of these factors can be 
made only if a knowledge of the motions produced by the 
controls is available. 

The angle of bank produced, for example, in 1 second 
by the total deflection of the ailerons might be taken as 
a measure of the effectiveness of the ailerons in produc- 
ing a sudden roll. If, on the other hand, the airplane 
is rolled slowly, a greater yawing motion being permitted 
to develop, the result may be to produce a rolling motion 
opposed to that generated by the ailerons and of such 
magnitude that the effect of the ailerons is completely 
nullified. The airplane may even roll against the ailer- 
ons. 


In actual flight the pilot might counteract this ad- 
verse rolling tendency by coordinating the rudder with the 
ailerons. In the present analysis, however, only the in- 
herent dynamic characteristics of the airplane are con- 
sidered. It would be of interest, then, to know the varia 
tion of the angle of bank with time after small aileron 
deflections are applied. Tor similar reasons, the effec- 
tiveness and the sensitivity of the rudder can be most 
carefully Judged if the angle of bank and the angle of yaw 
produced in a definite time interval by a small rudder de- 
fleotion are known. Information of this type is particu- 
larly valuable in estimating the effect of various combine 
tions of vertical fin area and wing dihedral on the dynam- 
ic behavior of an airplane. 
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and by the rolling and the yawing velocities) may be writ- 
ten: 



In order to solve for any of the variables, it is 
necessary to integrate this system of linear simultaneous 
equations. Tot reasons that will be apparent later, it 
is convenient to solve the system of equations (l) for 
separate unit magnitudes of the control disturbance terms, 
and 6Ng; that is, one set of solutions is obtained 
by letting 6Lg = 1 and 6Ng = o, and another set of 

„ * o uy letting 6Lg * 0 and 6Hg * 1. 

The unit disturbances are assumed to be instantly applied 
at zero time and to remain constant thereafter. In order 
to distinguish the separate solutions, the subscript L 
will be applied to the solutions obtained when 6Lg * 1 

and 6Hg * 0; whereas the subscript N will be applied 

to solutions obtained by letting 5Lg = 0 and 6ffg * 1. 

Thus, p L represents the rolling velocity resulting from 

the application of a pure rolling disturbance of unit 
magnitude. 

If the symbol D is substituted for d/dt and if 
6Lg =* 1 and 6Ng = 0, equations (l) may be rewritten in 
the fo Hawing form: 
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p (D-Lp) - rL r - PLp ** 1 
p (-Vp) + r(D-V r ) - fcVft = 0 


( 2 ) 


,(-£-) +I >, + D( B -It)p = 0 


It can he shown that the solutions of equation (2) 
are of the forint 


PL 


PL 0 + PL l * Xlt 4 p L a 0Xat * p L 3 ® X3t + p L 4 * A4 


X. t 


(3) 


where 


p- resultant rolling Telocity due to unit 
L rolling disturbance, that is, 6Lg = 1 


Pt --Pt 
L O L 4 


constants that depend only on values of 
stability derivatives and on type of 
disturbance involved 


roots of stability equation ?(D) * 0, 


where 



The constants in equation (3), together with the eoa- 
stents appearing in the expressions for all the remainiai 
components of motion (including the solutions for the unit 
yawing disturbance 6Vg = 1), have been evaluated in tsrs* 
of the stability derivatives and the roots of the etablUtf 
equation. The solutions of equation (2) are tabulated is 
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the appendix in a form suitable for computation. The 
stability equation is also discussed and an alternative 
form of equation (3) is given for use when the solution 
of 1(D) = 0 includes conjugate complex roots of the form 

a ± lb. 

The expressions presented in the appendix were eval- 
uated by applying the operational mathematics of Heaviside 
to equation (2). For the theory of operational methods 
the reader is referred to a standard text on the subject, 
for example, reference 1. Specific applications of the 
Heaviside treatment to other problems in airplane dynam- 
ics may be found in reference 2, 


After the complete unit solutions have been obtained, 
time histories of the motion caused by the unit disturb- 
ances can be plotted with very little additional calcula- 
tion. Because of the linearity of the equations of motion. 


the unit solutions may be compounded in any arbitrary man- 


ner. If, for example. 


8 a Ii£ represents the rolling ac- 


celeration created by the applied aileron rolling moment 
and 8 a Hg represents the accompanying yawing accelera- 


tion, then, at time t. 


0t = 6 a L 5. (O* * 5 a*8. (y*) t 


where 


resultant angle of bank after t seconds 

\^ l ') aagl® of bank after t seconds due to a 
^ unit rolling disturbance 

angle of bank after t seconds due to a 
' unit yawing disturbance 


er t seconds due to a 
tur bance 


As in the case of the unit disturbance terms, the 
actual disturbances, 8 a I*g and 6 a Yg * ar * as® 11 ®** 


be suddenly applied at sero time and to remain constant 
thereof ter. 
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App licat long*- The practicability of using the solu- 
tions to the lateral equations of motion is best illus- 
trated by concrete applications. 

The airplane under consideration, as mentioned in the 
introduction, exhibited certain undesirable flying charac- 
teristics in flight tests: The dihedral effect was un- 

desirably low at slow speeds with the flaps down (fig* 3(b) 
indicates satisfactory dihedral characteristics with the 
flaps up over the range considered), and it was impossible 
to raise a wing by use of the rudder alone. In an effort 
to improve the lateral flying qualities of the airplane, 
it was proposed to increase the wing dihedral. Because 
the lateral characteristics of an airplane depend not only 
on the absolute amount of dihedral but also on the relative 
amount of weathercock stability present, it was considered 
necessary to increase the vertical tail area as well as 
the wing dihedral. The ailerons were unchanged, but the 
rudder was so modified as to improve its hinge-moment 
characteristics. These modifications are Bhown on figure 
1 . 

Wind-tunnel tests were made with a model equipped for 
power-on tests with and without the proposed modifications. 
The data from the comparative tests indicated that con- 
siderable improvement should result from the incorporation 
of the modifications on the full-scale airplane. Nor all 
flap and power conditions the dihedral effect d C d\|/ re- 
mained positive, the index of weathercock stability dC n /d^ 

remained negative, and the static characteristics of the 
rudder appeared satisfactory. Wind-tunnel tests of a 
model, however, provide no direct information pertaining 
to the dynamic flight characteristics of the full-scale 
airplane. Instances have occurred in which the incorpora- 
tion of similar modifications on a full-scale airplane haa 
affected the control characteristics of the airplane in an 
adverse manner, particularly at high speeds, in spite of 
the favorable static characteristics indicated from wind- 
tunnel tests. Accordingly, it was decided to investigate 
the response of the airplane, with and without the modifi- 
cations to the aileron and the rudder controls, by evaluat- 
ing the solutions to the equations of motion. 

Aileron control.- In the particular problem considered 
here, undesirable aileron control characteristics are like- 
ly to be manifested in the form of aileron "heaviness" or 
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stiffness at high speeds. The subsequent analysis, how- 
ever, is perfectly general and is in no way limited to 
this high-speed condition. This particular condition is 
treated, merely as an example, to indicate the general 
method of procedure. Because the ailerons are identical 
on both the original and the modified airplanes, aileron 
heaviness can be physically interpreted as an increase in 
stick force resulting from the increased aileron deflec- 
tion necessary to reproduce a given rolling maneuver with 
the modified airplane. This interpretation of aileron 
.^saviness suggested the following method of analysis: 

(a) On the assumption that the original airplane 
(airplane A) was flying in steady high-speed flight, the 
angle of bank generated in 5 seconds by a small aileron 
deflection 8 a was computed. 

(b) Then, the aileron deflection 5 a + A5 a necessary 

to reproduce the identical maneuver with the modified air- 
plane (airplane B) was calculated. 

The time for the maneuver (5 sec in this case) Is 
somewhat arbitrary. It should be of sufficient duration, 
however, to permit full development of the secondary roll- 
ing effects introduced by the induced yawing motion. 

Tll0 Q i t lldO Of t. V\ • < nj'Tflmanf A St 4 — i i — a 

_ _ * u » v*. x x ov v 

measure of the additional stick force required to perform 
the maneuver (because the ailerons are identical) and may 
be used as an index of aileron heaviness. If A5 is 

large and positive, the aileron stick forces on airplane 
B may be too large to be acceptable and modification of 
the aileron itself may prove desirable. 

The stability derivatives of airplanes A and B 
were evaluated from the power-on wind-tunnel tests of the 
two models and from data in references 3 and 4. These de- 
rivatives are tabulated in table I, together with other 
information necessary to evaluate the solutions to the 
equation* of motion. 

In accordance with equation (4), the angle of bank 
assumed by airplane A in t seconds after the applica- 
tion of an aileron deflection 6. is 

& 
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" 8 » 1 6. (*l) t * 6.»6 a (*v\ 

In order to evaluate the unit solutions and 

it is necessary to obtain the roots of the stability equa- 
tion F(D) * 0 . If the appropriate values from table I 
are substituted into the expression for F(D) given in ths 
appendix. 

1(D) » D 4 + 20.4555D 3 + 52.7884D 8 + 347.8242D ♦ 5.43760 = 0 

With the use of the procedure outlined in the appendix 
for solving quartie equations, the roots of this equation 
were determined to be: 


= -0.01567 
X a « -18.6230 
X 3 * -0.908424 + 4.21991 
X 4 = -0.908424 - 4.21991 


Substitution of the appropriate roots and derivatives 
in the expression for the angle of bank given in the 

appendix, yields 


ft _ . ... _ _ . — o.oiS67t —is. east 

9^ * 3.3471 - 3.3497* + 0. 002876a 

+ f- O. 46778 + 0.1 6 310 93l\ (-0.90®* + 4.3i99i)t 
\ 1741.79 - 2270.451 /* 

+ f' -O « 46778 - 0 . 1 6310 931^ (-0*9084 - 4.3i9*i)t 
\ 1741.79 + 22 767451 /* 


If the last two terms are combined in accordance with 
the transformation formula in the appendix, the final ex- 
prootlon for 0^ can bo written: 

u „ _ ^ “ 0 . o 1 BS7 1 — east 

* 3.3471 - 3. 3497e + 0.002876e 

- 0 .00034623*"*° * 9064 ^ c0 9 4,2199 (t 4 0.1354) (5) 
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S imllarly , 

4^1 « 11.5982 - 11.62275e"°*° 1B67t + -0 .00038854e~ X 8 * 623 t 
+ 0 • 0 43548 e~ ° ’ 9 084 " C0B 4,2199 (t + 0.23305) ' ( 6 ) 

Practical solutions are moat conveniently obtained 
by graphical addition and subtraction of the component 
parts of the motion. It can be seen that a subtraction 
involves the small difference of relatively large quanti- 
ties. For this reason It is necessary to retain as many 
absolute figures as possible in the evaluation of the 
roots of the stability equation and in the evaluation of 
the individual components of the motion, 

Equations (5) and ( 6 ) are plotted in figure 5 (a) 
for airplane A. The corresponding expressions for air- 
plane 3 are plotted in figure 5 (b) f 

From figure 5 (a) the unit solutions for and 

after 5 seconds are 

$1 = 0.25 radian 

3 3 6 C 

H 0 30c = °- 66 radlan 


For a 1° aileron deflection from table 
6 a £g *» l f 54 per second 3 


and 


8 a » 6a * 0 


Hence, for airplane A, 

($) 5 #ee = (1,54) (0,25) m 0,386 radian =* 33, l 9 

The total aileron deflection necessary to bank air- 
plane B 32,1® i* 5 seconds is given by 
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In view of the fact that the ailerons on the two air- 
planes possess the same hinge-moment characteristics, it 
can he concluded that the aileron stick forces developed 
during maneuvers will at least he no greater (or heavier) 
on the modified airplane than on the original airplane. 

Rudder cont rol .- The sensitivity of the rudder con- 
trol, particularly at high speeds, may he judged hy the 
angle of hank and the angle of yaw generated, after a def- 
inite time interval, hy the yawing moment impressed hy 
the rudder. Here again the time must he of sufficient 
duration to permit the interaction of secondary rolling 
and yawing effects. 

Airplane A was assumed to he flying level at high 
speed and to he suddenly subjected to a yawing moment 
impressed by a given rudder deflection. The angle of hank 
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<j and the angle of yaw \|/ , generated after 5 seconds, 
were calculated. Under similar circumstances, the angle 
of hank 0 + A0 and the angle of yaw \J/ + A^ developed 
*y airplane B were similarly calculated. The magnitudes 
of A0 and A^/ were taken as a measure of the relative 
sensitivity of airplanes A and B to their respective 
rudder controls. The computations follow. 

From figure 5 and table I, 

Ai rplane A A irplane B 


W. ... ” dl,n 

CVJ 

. 

0 

1 
i 
i 

0.26 

(*). ... ” dl * n 

■ - - - .86 

.98 

6 r lg . per sec 2 - - - - - 

j. 

• - - - .308 

.410 

8 r Ng^, per sec 2 - - - - - 

■ - - - .549 

.652 


If a form similar to that of equation (4) is used. 


d. 

r t; 


/k t. . \ f d ^ ft vr \ f \ 

V w r-5 r J V L h VT»A 


and the angle of bank of airplane A, 5 seconds after the 
sudden application of a 1° rudder deflection, is 


8ec * (-0.308) (0.25) + (0.549) (0.86) 
* 0.395 radian = 23.6° 

For airplane B, 

#ec = (-0.410) (0.25) - (0.652) (0.98) 
= 0.532 radian = 30.5° 


and 


A0 = 8° 
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The angles of yaw generated hy the unit disturbances 
were obtained in a similar manner: thus, 

(+) * s r L 8 r (*l), + 5r * 6 r ('*'»). 

t t • 

In the maneuver investigated, the unit solutions after 
5 seconds were as follows: 


(O 

co 


radian 


s sec 


s sec 


radian 


Airplan e A Airp lane JB 
- 0.04 0.04 


.30 


.26 


Tor airplane A, 

(^) = (-0.308)(0.04) + (0 . 549) (0 . 30) 

5 8 0 C 

= 0.153 radian = 8.8° 
and, for airplane B, 

(dr) = (-0.410) (0.04) + (0.652) (0.26) 

s sec 

= 0.154 radian = 8.8° 


or 

A\| f 33 0° 


It appears that, for equal rudder deflect ione, the 
modified airplane will tend to generate more bank than 
the original airplane. The increase in tne wing dihedral 
and the modified fin and rudder may result, then, in a 
rudder control that will be slightly more sensitive at 
high speeds than the rudder control on the original air- 
plane. 


Hinge-moment measurements are, of course, necessary 
to determine whether the resultant rudder-pedal forces on 
the modified airplane will be greater or less than those 
on the original airplane. The hinge-moment characteristics 
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ean be compared on the basis of the quantity 



Tor airplane A, B * -0.297; for airplane B, R 
= -0.150. The resultant pedal forces, on the other hand* 
depend not only on the hinge-moment charact erist ice of 
the rudder hut also on its floating characteristics and 
on the particular deflection required to perform a stip- 
ulated maneuver. The resultant hinge moment per degree 
of rudder deflection may he expressed as follows: 


* 

6r 




The quantity 


'36 r \ 

.■§*70* = 

r 


0 


is theoretically a con- 


stant for any given tail arrangement hut actually it is 
very critical to interference effects at the tail and ^ 
fluctuates considerably. Tor small angles of yaw 5 j , 



Ch 


r 


was practically sero for both airplanes A 


and B; and the quantity R may,, therefore, he taken as 
a measure of the pedal forces for a given rudder deflec- 
tion. 


TBS LOSQXTUDXSaL MOTION OT THB AXRPLANI 


then the longitudinal stability of an airplane is 
discussed, the characteristic usually referred to is the 
"static" longitudinal stability. If static longitudinal 
stability exists, the dynamic stability characteristics 
are of minor importance (reference 5). 
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The usual index for static longitudinal stability is 
the rate of change of pitching-moment coefficient with 
lift coefficient dC m /dCi, or some quantity proportional 

to it. In flight the most convenient method of evaluating 
the amount of static stability present is to measure the 
elevator angle required to trim the airplane at various 
speeds; the slope, dS e /dV or d6 0 /da, is an index of 

the degree of static stability possessed by the airplane. 
The significance of the ratio d6 e /da. methods of eval- 
uating it for power-off and windmilling conditions, and 
suggested design values are discussed in reference 6. The 
necessity for power-on wind-tunnel tests for securing the 
effect of power has already been discussed (fig. 2). 

If tunnel data are available from which the floating 
angle of the elevator may be ealuclated for any lift coef- 
ficient, the elevator stick force required for trim at any 
lift coefficient can be calculated from the following 
formula: 


where 


P 






( 7 ) 


P stick force for trim, pounds 
x linear travel of top of control column, feet 

A6 0 difference, in degrees, between elevator 
angle required for trim at lift coef- 
ficient under consideration and free- 
floating angle of elevator 

The term A5 e ° may fluctuate considerably with power. 


The slope of the curve relating the variation in stick 

v 4 C v !v th forward «P e « d d P/dV depends on the speed at 
wnich the airplane is trimmed for zero stick force and, 
consequently, depends on the initial setting of the trim- 
ming tab. In the comparison of the change in the elope of 
the stick-force curve resulting from modifications to the 
elevator, care should therefore be taken to orient the 
elevator trimming-tab settings so that zero stick force 
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always occurs at the same speed. Otherwise, a superficial 
examination of the curves of stick-force variation with 
air speed may lead to incorrect conclusions concerning the 
effectiveness of the tail surfaces or the. stability of the 
airplane. 

In addition to providing a means of trimming tha air- 
plane in steady flight, the elevator must be capable of 
changing the airplane flight path. The rate at which this 
change is accomplished in a quick pull-up can, in a sense, 
be Interpreted as a measure of the effectiveness of the 
elevator in maneuvers. 

It is convenient to take the rate of change of the 
maximum normal acceleration per unit of elevator deflec- 
tion as an index of elevator effectiveness in maneuvers. 
Care must be exercised, however, in interpreting this 
index. Although it is essential in a pursuit airplane to 
design an elevator sufficiently powerful to maneuver the 
airplane to the maximum lift coefficient of the wing, it 
has been found very undesirable if this condition is ful- 
filled with a minimum amount of elevator deflection. As 
discussed in reference 6, satisfactory static stability 
characteristics require the quantity dS/da to have a 
value around 0.5. In airplanes that required consider-' 
ably less stick travel to trim the airplane over the 
angle-of-attack range inadvertent stalling has frequently 
occurred in accelerated maneuvers. The optimum value for 
the rate of change of normal acceleration per unit of 
elevator deflection is therefore conditioned by the re- 
quirements of satisfactory static longitudinal stability. 


In the following section a simplified formula is 
developed that permits the rapid estimation of the normal 
accelerations developed in abrupt pull-up maneuvers. 


D evelopmen t of a simpli fi ed form ula for normal ac- 
celerations produced in abrupt -pul l -ups .- To a first 


order of approximation, the 
plane of symmetry involving 
written as follows* 

Sr 3 u^ + w*w + «*, - 

= * WZ* + q (u 0 

It* = uM u + wM w + + 


equations of motion in the 
a disturbance in pitch may be 


ge 

(a) 

+ 2q ) - sB 0 e + Se z 5 e 

(*> 

6e«5 

(c) 
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Also 


d0 _ 
dt ” 


q 


where 9 0 is the initial angle the X axis makes with 
the horlsontal axis. 

The axes in which equations (8) are expressed are 
similar to those defined in the section on lateral motion 
except that the X axis is inclined at an angle 0 to the 
relative wind. 

Xn general, this system of equations must he solved 
by methods analogous to those used in obtaining the solu- 
tions to the lateral equations of motion. The normal ac- 
celeration is given by the expression dw/dt - qU 0 . 

These equations have been solved for the airplane, 
the characteristics of which (evaluated, so far as pos- 
sible, from wind-tunnel testB on a powered model) are 
presented in table II. The curve of the normal accelera- 
tion resulting from a 1° upward elevator deflection is 
plotted in figure 6 together with other components Of the 
motion. 

The variation of the components of motion shown in 
figure 6 is typical of an abrupt pull-up maneuver at high 
speed. The following facts are apparent from this figure: 

1. During the time the airplane takes to attain 
maximum acceleration, the velocity V remains sensibly 
constant and is equal to U 0 , the equilibrium velocity. 

Accordingly, u *= du = 0 , and equation (8a) can be ne- 
glected. 

2. In the vicinity of the maximum acceleration, the 
time rate of change of pitching velocity dq/dt is ap- 
proximately sero, that is, q is approximately constant. 

3. In the vicinity of the maximum acceleration, the 
acceleration component dw/dt is almost zero and the ac- 
celeration thereafter is given almost entirely by the 
product qU 0 . 
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If Z 5 


and g0 o (all of which are small) are 


^ "q • o 

neglected and the substitutions CtZ a = wZ, and »*w 
= aM a are made, equation (8), for an abrupt pull-up ma- 
neuver, can be simplified to the following form: 

qU 0 + aZ a =0 ^ 


qM q + OM a 


-M, 


where 


S e*6 


e 


The eolution for qU A is 


where 


W 0 


"o 


/ “ 0 Z a \ 

^o“a- 2 a«q) 


Mo = 


Z a = 
M a = 


££m\ . 1 P5^ 8 c 

36 J 5 ' 2 hoc/ 

(Ml) i p sv 3 

\da/ 2 m 

dC m \ 1 P ST" c 

<da J 2 


“q 



1 PSVI 3 

2 mky 2 S 


(9) 


In the expression for Mq , the slipstream factor T 

normally has a value between 1 and '*2 5 for high-speed 
flight. The tail efficiency factor r) t is always less 

than 1 and is generally about 0.9. The value of the pro* 
uct *Dt 18 therefore always about unity. 
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If the airplane is assumed to "be flying level before 
the pull-up maneuver, mg = C L ipSV 2 = CjipSUQ 2 . If the 

C C 

euhst itut ions U 0 = V , p, = r-® — , and = 1 are made, 

2 PSI 

equation (9) can he reduced to the following simplified 
form? 

Q^o _ _ S«t dOm/ dSg. 

5 a C L u dC m _ l 7 * Cl 

do i c \da 

where qU 0 /6 # represents the change in normal accelera- 
tion per unit of elevator deflection. 

Formula (10) gives a value of the normal acceleration 
produced during an abrupt pull-up maneuver that closely 
approximates the maximum acceleration which would be ob- 
tained by solving the more cumbersome equations of motion. 
It is the acceleration that will be obtained if the ele- 
vator is instantly deflected to its final position and 
held in that position until the maximum acceleration is 
reached. The formula gives less accurate results when 

dCm/dCj, is small, that 'is, - — S. <0.01. Several cal- 

dc L 

culations made for smaller values of dC^/dC^ gave re- 
sults, however, that were in error by less than 7 percent. 
Greater accuracy is to be expected at small values of 

on account of the approximations made involving 0. At 
high values of Cjj , these approximations introduce greater 
errors. 

M^hfUTorabi l ity and stabil ity .- The manner in which 
the normal acceleration produced ir a pull-up maneuver is 
affected by d^/dC^ and p, is shown in figure 7. An 

examination of this figure reveals that, for airplanes 
with high wing loadings (large p,) and low static stability 
( small dOm/dCL) , the normal acceleration produced per 
degree of elevator is affected by a small change in the 
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static stability. In accelerated maneuvers this effect, 
as far as the pilot is concerned, will manifest itself in 
the form of increased stick forces if it is assumed that 
the pilot wants to produce a given acceleration regard- 
less of the degree of static stability present in the air- 
plane. 

The index of static stability dC m /dC L for an air- 
plane with the characteristics given in table II is as- 
sumed to be -0.022. The curve of the normal acceleration 
produced in a pull-up from level flight at 448 feet per 
second is given in figure 6. The maximum change in the 
normal acceleration attained with an upward elevator de- 
flection of 1 is about 84.5 feet per second. 

From formula (10) 

lEo. = ... 3 4 « 4 ( -0. 021 > \ _ _ c 

5 q ° 6 0.10 \( 34 . 4) (-0.022) + (-2.01)^ = 2 * 6g P er 

degree and qU Q = 83.7 feet per second 2 , an error of less 
than 1 percent. 

If the static stability of the airplane is improved 
by moving the center of gravity forward 0.078c, the ex- 
pression for the normal acceleration becomes 


qU s — g S. 4 * 4 —0.021 e 

0 0.10 \( 34. 4) (-0.10) + (-2.07)/ 


= (- 1 .37g) 5 # ° 

If it is desired to load the airplane to its former 
additional load factor of 2.6g, the elevator movement re- 
quired with the higher dC^dC^ would be 


5 


e 


= -!.9° 

- 1 . 37g 


Accordingly, the pilot would have to exert nearly 
twice as much stick force to execute the identical pull- 
up maneuver because of the increased static stability of 
the airplane. The desirability of modifying the elevator 
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in order to change its hinge-moment characteristics thus 
partly depends on the relationship between the static 
stability desired and the magnitude of the stick forces 
acceptable in accelerated maneuvers. 


CONCLUSIONS 


This paper is intended to illustrate primarily how 
power affects the characteristics of a model tested in the 
wind tunnel and how wind-tunnel data may be used to esti- 
mate flying qualities not directly discernible from wind- 
unnel tests. The analyses presented in this paper permit 
the following conclusions to be drawn concerning the meth- 
ods employed: 

(1) In the prediction of the flight characteristics 
of an airplane operating with power, considerable error 
may be introduced if wind-tunnel data from tests of a 
model not equipped with an operating propeller are used. 

(2) In the analyses of wind-tunnel rolling-moment 

a a, care should be taken, in the determination of the 
ihedral effect, to allow for the contribution of the un- 
balanced pitching moment to the slope of the rolling- 
moment curve. 3 

(3) The evaluation of the equations of motion permit 
estimates to be made of the relative effectiveness of the 
ailerons and of the rudder controls, particularly when 
changes in wing dihedral are involved as these effects are 
not readily diacernible from static tests of the controls 
alone. The methods used, have been found reasonably ac- 

wllen "ind-tunnel data are available, and are not 
difficult to employ. 

^ (4) The rate of change of normal acceleration per 

unit of elevator deflection affords a convenient correla- 
tlon between maneuverability, stability, and elevator- 
stick forces. The approximate formula developed for cal- 
culating the normal accelerations produced during an 
abrupt pull-up maneuver is simple to evaluate and has been 
found to yield reasonably accurate results, particularly 
at high speeds. 


Langley Memorial Aeronautical Laboratory, 

National advisory Committee for Aeronautics, 
Langley Field, Va. , July 22, 1941. 
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APPENDIX 
List of Symbols 

velocity along axis in steady flight 

velocity along flight path 

velocity along X axis 

sideslipping component of velocity 

velocity along Z axis 

angnlar velocity in roll 

angular velocity in pitch 

angular velocity in yaw 

angle of hank 

angle of yaw 

angle of attack 

angle of Bideslip (®v/ T J 0 ) 

angle X axis makes with horizontal 

control setting, with appropriate subscript 
(6° indicates that values are in degrees 
and not in radians) 

components of force along X, Y, and Z axes, 
respect ively 

rolling moment about X axis 
pitching moment about Y axis 
yawing moment about Z axis 
hinge moment 
wing span 

chord (of wing, unless otherwise subscripted) 
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S area (of wing, unless otherwise subscripted) 

l tall length (distance from center of gravity 
to tail post) 


m mass 

u relative density factor — -i£ 

l/2pSl 

mkj 2 moment of inertia about X axis 

mky a moment of inertia about Y axis 

mkg 2 moment of inertia about Z axis 

t t lme 
D d/dt 



density 

gravity (32.2 ft/sec 3 ) 


ro 1 ling- moment coefficient 


( 

Vi/SpV'Sb 


) 


pitching-moment coefficient 



C n yawing-moment coefficient 
Cfc hinge-moment coefficient 
Cy lateral-force coefficient 



0^ lift coefficient 

T) t tail efficiency factor 


f empirical slipstream factor to account for con- 
tribution of fuselage, wing, and propeller to 
damping in pitch 
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lat eral- stability derivatives in terms of unit 

mass or moment of inertia of airplane. 3Por 
example. 


Y 6 = 
r m 

L_ = ^./ 5 ,P 

"V 

u = 

T 



longitudinal- stability derivatives in terms of 
unit mass or moment of inertia of airplane. 
For example. 



M u = M/du 

mky 3 
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Subscripts are defined: 

L solutions obtained when SLg = 1 and 6Ng = 0 

N solutions obtained when 61*5 = 0 and 6Ng = 1 

a aileron 
e elevator 
r rudder 
t time 

o Initial value 

Primed symbols refer to horisontal tail. 

SOLUTIONS OP TEX LATERAL EQUATIONS OP MOTION WHIN 
THE ROLLING DISTURBANCE IS UNITY AND THE YAWING 
DISTURBANCE IS ZERO (5Lg = 1; 6Ng = 0) 


The stability equation P(D) =0 is obtained from 
equation (2) of the text by expanding the third-order de- 
terminant formed by the coefficients of the variables. 
Thus, 
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-10 

* ~(h + »r + J&V 

+ 57 (*i- + 1 p) + ^rlp * Vp + Sp 0 2 

+ - ^ <S l P * l rSp> S K P + ^ B p -£»,»» 

+ 5 - (lf* r - lr B |s) = 0 



Let A x » A 2 , A 3 , and A 4 he the four roots of this 
equation J(D) = 0 and form the following products: 


p ■ \ X 3 \ \ = (l 6 JT r - L r Rp) 

« » X l ( X l - * a ) (*i— X a ) (x, - X 4 ) 

8 - X E (x a - \) (x a - X 5 ) (x a - X 4 ) 

8 - *3 (*» - \) (8. - 8.) (x, - X 4 ) 


I * X 4 (x 4 - X,) (x 4 - X a ) (x 4 - X 5 ) 

«ra a «!?*/°* nti0n8 for any of th ® var ia^les can be ex- 

and the °5 *?* products of »he roots just formed 

d the stability derivatives. By use of the auxiliary 
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relations 7 T = * and r~ = r solutions can be written 
at at 

directly for both 0 and \j/. As it is usually desirable 
in this work to determine 0 and rather than p and 

r, it is convenient to solve for 0 and V directly and 
then to differentiate each of these solutions, respective- 
ly, to obtain p and r if these variables are wanted. 
The solutions follow: 


Angle of Bank 0^ 
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Angle of Sideslip 


h 


where 










4 


X 3 t X 3 t 

+ p l . + p L • 

Uf J 


- g -N r 
u o P 


_ g 

(i>- ^ 

)-_ x £e 

0 ° 4 

1 

k 

(x 2 - u r/ 

) ' 

U 0 

R 


_ 6 

0= - *r' 

) - S»p 

U 0 

s 


_ _g_ 

(*4 - Bry 

Lr-iia 


U 0 T 


+ 


Angle of Yaw \j/ L 


X„t, X,t 

^1 3 'J't, + 6 + e ° + \|f- e + 

O 1 2 "3 


where 
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* 


1 


o 


_g 

P 




l 






3_ 




Vp 




VP 

c o 


Tawing Velocity r^ 


r L = i (+ L ) 


The yawing velocity is most easily obtained by direct 
different iat ion,. 


Rolling Velocity 



The rolling velocity ie most easily obtained by di- 
rect differentiation. 
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SOLUTIONS OF THE LATEEaL EQUATIONS OE MOTION WHEN 

TH1 yawing distubbance is unity ahd the bolling 

DISTUEBANCE IS ZEEO (5N g = 1; 6L g = 0) 

Angle of Bank 




X.t 


^N * + ® 1 + $N _ e 2 + ^TB ® + ^N * 4 

o l "e 3 "4 


where 


l£ 

_ ~H “ L r U 0 
P 


A 


K - l r S + \ 


N. 


“Lc - 


5 N 


*ai 

L r + Vr 


E 


[-If - L r It) + X 3 L r 


’F. 


J 




^ * Nlr 
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Angle of Sideslip 


X.t X 2 t 


V 


f>* = h 0 + *>!,• 1 + ^ 2 e + P »3* + 


V 


where 


gXr 

b - J:°- 
** 0 " f 


ibl + XL- X £ 
U ft IP 1 


3 17 


~ + X L»j - Xj 
U 0 2 P 3 


^3 = " 


gL 

-- + X 3 Lp - 


** - 


fr ♦ V P ' V 


Angle of Taw 


B 


■ **. + V M + + + 


X„t 


x,t 


X.t 


HACA Technical Note No. 828 


35 


where 




g 

>L & u<T 


* 


x 2 Y p Yp 


•- L 


J 

p ur 


i. 


= 


->.„ 3 !& + X„L. Is.ui. 

_ 1 _ U 0 __ ®_P U© ^_Uo 

R 


V 


, = H , T P « 

3 Uo + 3 L p u 0 " L p Uo 





*4*p 

~T 



Yawing Velocity rjj 


r R 


d 

dt 



The yawing velocity is most easily obtained by direct 
dif f erent iat ion. 


Rolling Velocity pjj 

p * = A W 

The rolling velocity is most easily obtained by 
direct differentiation. 
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DI80USSI0S OF THE K 4 UATIOU F(D) = 0 


For most airplanes the solution of F(D) = 0 will 
yield both real and imaginary roots. In the solution of 
the equation, the real roots can first he isolated 
(Horner* s method) and extracted from the equation. The 
Imaginary roots can then he found hy solving the result- 
ing quadratic. Because imaginary roots always occur in 
pairs# two of the roots will he of the form a ± ih, 
where a and h are constants. 

The components of the motion containing the imaginary 
roots can he combined conveniently into a single term in- 
volving only real numbers. In the case of the solu- 

tion is 


<*i ■ *i„ + *i • 


X^t . Xgt j ^3^ . J 

+ h • + K • + 


If X and X. are the conjugate imaginary roots. 

3 ' 


x,t x 4 t 

i + 0 . e will also he imaginary. 


let X, * a + h; then 


will he imaginary and 


can be reduced to the form ?L_£_i.£.. Further# 

0 + iH 


S_t_iF iH - SO. + _FB .- i (BB._r.JgOl = L + 

& + iE 0 - iH a 2 + H 3 1 


and 

X t X t 

h* * - *!«• 


t A t o / > S i 

3 + dm m = — e at cos h (t + y) 


where 


7 = g tan 
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rig*. 2,3 
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figure 2.- Variation of elevator angle required for trln with Indicated airspeed . 
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